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We study a family of local boundary conditions for the Dirac problem corresponding 
to the continuum limit of graphene, both for nanoribbons and nanodots. We show that, 
among the members of such family, MIT bag boundary conditions are the ones which are 
in closest agreement with available experiments. For nanotubes of arbitrary chirality sat- 
isfying these last boundary conditions, we evaluate the Casimir energy via zeta function 
regularization, in such a way that the limit of nanoribbons is clearly determined. 
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1. Introduction 

It would be redundant to start this paper with a detailed description of the wonder- 
ful properties of graphene. In this respect, the interested reader is referred to Ref.[I] 
As predicted theoreticallj^Sl twenty years before its production in a laboratory, 4 
electron transport in graphene is d escr ibed by a massless Dirac equation, which leads 
to distinctive electronic properties Indeed, many such properties were studied ex- 
perimentally, starting with the determination of the quantum Hall effect,^ and found 
to agree with the predictions of a "relativistic" and massless Dirac field theoryP 



•Member of CONICET 
t Member of CONICET 



December 23, 2011 1:22 WSPC/INSTRUCTION FILE qfextbeneventano 



2 C.G. Beneventano and E.M. Santangelo 

The aforementioned properties allow to envisage many possible applications. 
However, a crucial point in achieving such goals as the construction of graphene- 
based transistors is the opening of a controllable band gap in an otherwise gapless 
material. The use of samples of finite size is a natural guess when trying to do sof^ 
In fact, several measureme nts of the electric conductivity in graphene devices do 
show the existence of a gapJHUIIIIIIHI 

Unlike the case of usual semiconductors, the confinement of charge carriers to a 
finite region cannot be modeled, in the continuous Dirac theory, by the condition 
that the fields vanish at the boundaries. Most theoretical approaches to this pr oblem 
presuppose an orientation dependence of the adequate boundary conditions) ^^ 
which is in contradiction with the experimental results.^ A nice general study of 
possible boundary conditions in the Dirac problem and of their symmetries in the 
case of nanotubes was presented in Ref. [15] 

The aim of this paper is twofold: in the first place, we will study a family 
of boundary conditions satisfying rather standard criteria, and compare the pre- 
dictions arising from taking different members of such family with experimental 
results. After having shown t hat , among these members, the ones corresponding to 
MIT bag boundary conditional^ are good can didates to reproduce the experimental 
outcomes, we will evaluate the Casimir energy ^^^ for nanotubes of arbitrary chi- 
rality satisfying such boundary conditions, and arrive to an expression which leads 
to a direct identification of the Casimir energy in the nanoribbon limit, and shows 
clearly the (exponentially decreasing) corrections due to finite length-to- width ratio. 

2. General setting. Edge states 

We will choose the orientation of the lattice as in Ref. [2j so that, by taking the two 
nonequivalent Dirac points as K± = (0, = ^ 3 ^ n ); we S e * the total Hamiltonian as a 
direct sum of 



H± = hv F (-ia 2 d x ± ioidy) , (1) 

where Vp — l^r is the Fermi velocity of graphene, with a = 0.14nm the distance 
between nearest neighbors and t = 2.7eV the nearest neighbor hopping energy. 

Such Hamiltonian corresponds to a free Dirac equation in 2+1 dimensions, where 
the gamma matrices are given, in each valley, by = 1173, 7± = 01, 7± = ±02- 
We will study the corresponding eigenvalue problems H± , $>±(x 7 y) — E± 1 & ±{x,y), 
when the domain of the differential operator is defined by a family of local boundary 
conditions which 

(i) Are separately imposed in each valley, 

(ii) Give a vanishing flux of current perpendicu lar to the boundary}^ 

(iii) Are defined through a self-adjoint projector fEUSl 



Note that the condition (ii) is as close to confinement as one can get in a Dirac 
theory and it leads to a self-adjoint Hamiltonian (thus, to real energies). 
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From now on, we will study the problem around K +1 leaving the discussion on 
how to combine boundary conditions in both valleys for later on. We will consider 
the boundary to be placed at a given x = xq value. Throughout our calculations, we 
will take vp = H = 1, and recover the right units when comparing our predictions 
with experimental results. 

We start our analysis with the case of the half plane, with its boundary at 
x = 0. For such geometry, it is easy to check that the current perpendicular to 
the boundary is proportional to ^\a2^+, while the current along the boundary is 
proportional to ^> + a\^! + . So, the most general local boundary conditions satisfying 
conditions (i) to (hi) above are given by (/ + o\ e _M<T2 )4 , +Jx=K = 0, which is a 
one-parameter family. In Ref. [T31 this family was associated to the existence of a 
staggered potential, different for A and B sites, in a total of 2N rows closest to a 
zigzag edge parallel to the y axis, as in our case. 

Note that a — 0, n correspond to the so-called MIT bag boun dary conditions}^ 
while a = ±^ are the conditions used to mimic a zigzag boundary!^ Here, we point 
out that, among the graphene community, MIT boundary conditions are usually 
called infinite-mass or Berry-Mondragon boundary conditions because these two 
authors studied them in the 2 + 1 context of quantum billiards, in Ref. [T9l 

Each member of this family of boundary conditions imposes a different con- 
straint on the density of tangential current at the boundary. In fact, one has 

ai^> + \ x=X[) — — cos (a^^^Z+lx^xg. In particular, zigzag boundary conditions 
enforce the tangential current to vanish at the boundary, while MIT ones equate it 
to the density of charge. 

In view of the translational invariance along the y direction we will propose, for 
each k y , ^ + (x,y) = e lky v ip + (x). In order to analyze the existence of edge states 
(or the lack thereof), it is convenient to perform a unitary transformation of the 
eigenfunctions, ip{x) — e~ l ~ a ' 2 ijj+(x). This leads us to the eigenvalue problem, 

\—i<Jid x — <J\k y cos a + a^ky sin a] tp(x) = Eip(x) 

(l + * 1 )#(x = 0)=0, (2) 

together with the normalizability condition when x — > oo. 

It is a simple exercise to show that, for all a ^ 0,n (i.e., all boundary condi- 
tions different from the MIT bag ones) there are, apart from bulk sates, some edge 
states, corresponding to E = fey cos a, with fc^sina > 0, which are eigenfunctions 
decreasing exponentially with x, thus concentrated at the boundary. 

The existence of these states is well known for zigzag boundary conditions (a = 
±f ), in which case they are zero energy modesP^ This shows that, in a compact 
region with a single smooth boundary, these two boundary problems will not satisfy 
the Lopatinsky-Shapiro condition (equivalently, they will not define a Fredholm 
operator).^" In the remaining cases, the edge states correspond to real, nonzero, 
energies. We will comment on this point when treating circular quantum dots, in 
section |U 
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3. Graphene nanoribbons 

For a comparison with experiments, we will first discuss the case of a graphene 
nanoribbon, which requires the imposition of a boundary condition at a second 
boundary, placed at x — W . _____ 

Experiments concerning nanoribbons^ _ __ snow a g a p which, moreover, is sym- 
metric around zero gate voltage. There are two ways of obtaining a symmetric spec- 
trum, i.e., choosing exactly the same projector to define the boundary condition at 
x = W or choosing, instead, the orthogonal one. It is easily shown that the first 
alternative allows for the existence of zero modes, no matter the value of a. They 
appear for all values of k y when a = ±-§ and for k y = in the remaining cases. So, 
we will limit our discussion to the second alternative, which is also consistent with 
the fact that the sign of the inward normal is opposite at both boundaries. 

Our boundary value problem will now be, 

[—i(J2d x — <J\k y cos a + a^ky sin a] ip(x) = Eip(x) 
{I + <t 1 )${x = Q) = Q 

{I-ai)i){x = W) =0. (3) 



In all cases one has E = ±y fc^ + ky. However, MIT bag boundary conditions 
are unique in that the spectrum is determined by the equation cos(k x W) = 0, 
which does not depend on k y , and they allow only real values of k x . Thus we have, 
in these two cases (a = 0, tt), 




E n = ±\ J' +kl n = 0,...,oo. (4) 



This leads to an energy gap = ^. We stress that the same dispersion formula 
was obtained from a microscopic model in Ref. 1211 

The remaining values of a, instead, lead to a spectrum determined by 

k x cos (k x W) — k y sin a sin (k x W) for E ^ ±k y , (5) 

and 

k v = T77-- J ioi-E = ±k y . (6) 

Note that both equations break the invariance under k y —ky. It is difficult 
to imagine why this invariance would be broken in a ribbon, which extends to 
— oo < y < oo, in the absence of electromagnetic fields. Such invariance could only be 
recovered by imposing exactly the same boundary conditions on the eigenfunctions 
around the other valley. 

(ft _|_ _. 

Moreover, for k y — one has k x = - — p^— , no matter the value of a. For all 
the remaining values of k y , at variance with the situation in the MIT case, the 
admissible values of k x are not equally spaced. 
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But, more important, Eq. ([5]) allows for imaginary as well as real values of k x . 
Calling k — ik x one has, for E ^ ±k y , 



k cosh (kW) — k v smas'mh(KW), for|fcj> — ; . (7) 

W sina 



When | sin a | = 1, i.e., zigzag boundary conditions, this equation allows for en- 
ergies arbitrarily close to zero when k — > oo. As a consequence, no gap exists in this 
case, which is a well known fact W For the remaining values of a, the eigenenergies 
coming from Eqs. ([5]) and ([7]) never tend to zero. The analysis of the minimal value 
of \E\ can be performed analytically. From such analysis one concludes that, for all 
a, the energy gap satisfies AE < 

The experiments on nanoribbon d^ l 1 1121 show a transport gap as a function of the 
gate voltage, when performed at low temperature and bias voltage. This eliminates 
zigzag boundary conditions as candidates to describe the physical situation. For the 
remaining values of a we have, recovering units, AE < hv ^ n = §7rf-^F- The equal 
sign holds only for MIT bag boundary conditions (a = 0, ir). 

As for the numerical value of the gap, Ref.[12]presents a study of several graphene 
nanoribbons of different widths, all of which show a gap in the gate voltage cor- 
responding to a one-particle energy gap fitted to A m = 30>eV^. This is roughly 
three times our result for MIT boundary conditions, i.e., AE = Yl.leV-^. Ref. IT01 
finds, for a sample of width W = 30nm, a value of the energy gap AE = A6meV, 
in much better agreement with our result for such case, i.e., AE = 57meV (note, 
that, in this case, our prediction is higher than the measured gap). Obviously, both 
experiments disagree. The origin of such discrepancy is not clear to us, since both 
use similar values of the capacitance for comparable samples. So, numerical values 
of the gap cannot be used, at present, to select among different non-zigzag values of 
a. However, Ref . [TU] shows equally spaced plateaux in the conductivity, which only 
MIT bag boundary conditions can explain. Thus, these are the best candidates for 
the boundary conditions to be imposed in the continuous model. 

Finally, as shown clearly by Ref. GO and at odds with most theoretical 
modelJISnS (which would impose different boundary conditions depending on the 
orientation of the boundary) the measured gap in the gate voltage doesn't depend 
on the orientation of the boundary. This will certainly be the case in our continuous 
model if MIT bag boundary conditions are written as (/ + fi)ip(x = Q,W) = 0, 
where n is the inward (or outward, which leads to the same result) pointing nor- 
mal vector corresponding to each boundary. Moreover, consistently imposing MIT 
boundary conditions around K_ leads to the same spectrum in both valleys. 



4. Nanodots 

In order to further compare the predictions of different members of our family of 
boundary conditions with the experiment, we will now treat the case of a circular 
graphene dot of radius R. To this end, we adopt polar coordinates. Taking the 
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gamma matrices for the theory around K+ as before, we are led to the boundary 
value problem (with vp = h= 1), 

r 

(i - Ye~ iayf> ) ip(r = R,6) = 

1 p(r,0)=^(r,9 + 2- ! r), (8) 

where r y r = o\ cos + o<i sin 9 and = <Ji cos 6 — cr\ sin 9. 

Solving this boundary value problem is a simple exercise (see, for instance, 
Ref. |2"2"]) . The first outcome is that zigzag boundary conditions (a = i^) allow 
for an infinite amount of zero modes, as expected from the facts that they do not 
satisfy the Lopatinski-Shapiro condition and that we are now treating a compact 
region with a smooth boundary. As for the remaining conditions in the family, none 
of them allows for zero modes. Since experiments on quantum dots also present a 
gap, we will concentrate on the remaining cases (cos a ^ 0), which give a spectrum 
determined by 

(1-sina) J n (\E\R)+s cos a Jn+i(\E\R) = 0, n = 0, oo 

(1- sin a) J n+1 (\E\R)~s cos a J n (\E\R) = 0, n = 0, ...,oo (9) 

where J n is the Bessel func tion of order n, and s is the sign of the energy. 

Now, the experimentPJJ shows a gap in a quantum dot which is symmetric around 
zero gate voltage. This, again, points to the MIT boundary conditions as good 
conditions to impose on the continuum model in order to reproduce the experimental 
results, since all the remaining values of a produce a spectral asymmetry, as can be 
easily seen from Eqs. ([9]). 



5. Casimir energy of graphene nanotubes and nanoribbons 

Now that we have justified the use of MIT bag boundary conditions when treating 
graphene devices in the continuum limit, particularly in the case of ribbons, it's 
interesting to calculate the vacuum expectation value of the energy, also known as 
Casimir energy, due to the finite size of a given sample. To achieve this goal, we 
will evaluate such magnitude, in the framework of the zeta function regularization, 
for a general graphene nanotube of arbitrary chirality (characterized by |), com- 
pactification length L and finite length W in the perpendicular direction, imposing 
MIT bag boundary conditions at the boundaries x = 0, W . The limit corresponding 
to graphene nanoribbons is — > oo. So, we will analytically extend the otherwise 
diverging sum of allowed energies in such a way that this limit is made patent, 
and the exponential corrections due to finite length effects are clearly displayed. A 
similar calculation was already presented in Ref. I23( in the particular case of an an- 
tiperiodic compactification (| = i). Here, we generalize the calculation to arbitrary 
< | < 5, which includes the interesting cases of conducting nanotubes (| = 0) 
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i). Due to the invariance of the spectrum under 



and semiconducting ones (| 
| -¥ | + 1 and | -¥ 1 - §, 
direction. 

The Casimir energy of the problem at hand, will be regularized according to 



this will cover all the possible twists in the compact 



E G 



9s 9v 



E E rJ+ E ^ 



E nJ >0 



E n ,i<0 



(10) 



where E n j are the modes in Eq. (jU, with k y = + to account for the different 
possible chiralities of the nanotube, and g s — 2 is the spin degeneration and g v = 2 
is the valley degeneration. 

More explicitly, for the Casimir energy per unit compactification length, we have 



L 



L 



OO OO 

E E 

l — — oo n—0 



. 1. 7T 



5 .2tt 



(ii) 



s=-l 



or, after Mellin transforming the previous expression, 

' o \ — S /*oo 00 00 r o ,i 

2n\ g s g v f ^i-iy- y- e -*{[(n+l)i^] 2 +[i+l] } 



L LT i 



I— — oo n—0 



(12) 



An analytical extension, convenient to isolate the term remaining in the nanoribbon 
limit, can be performed by making use of the inversion formula^ 



-i(/+c) 2 



1 oo 

-) 3 E e -^- 2 ^ c . 



(13) 



When doing so, we get 



Ec 
L 



OO OO 



z=i ra=0 



J2 / dti^e"'^^] .(14) 



Now, after performing the integral in the first term, and writing the second one 
as a Hurwitz zeta function (Ch), we obtain 



Ec 
L 



2ir \ s g s gv^ 2 



oo oo 



L J LT\ 



! 4 EE cos (^) 

1=1 n=0 



7T/ 



,( n + 2) 2W, 

s- 1\ L x _ 



2 /W™-^' 



.(15) 



s=-l 
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Finally, after relating the Hurwitz zeta function to the corresponding Riemann 
one (Cr), and using the reflection formula for this last pS we obtain the final ex- 
pression for the Casimir energy per unit compactification length, i.e., 



Ec 



2 cos(7r/d) 



LW 



irlL 



Z=l,ra=0 



W 



3g s g v 
32ttW 2 



Cfl(3) ■ (16) 



Not unexpectedly, the nanoribbon limit (^ — > oo) of this Casimir energy, given 



b y ^ - "IP 1 

Indeed, it could have been obtained, by considering a continuous k y in Eq. (|11[) . as 



- ^ Cfi(3)i turns out to be independent of the chirality of the nanotube. 



E c 
~~L 



9s 9v 
2tt 



dky ^ 

n=0 



(n 



-)- 



(17) 



s=-l 



The numerical value of this Casimir energy for nanoribbons is, after recovering 
physical units, 

Ec 
L 



._^ C h(3) = -0.1435-^. 



-0.2870§f. 



which corresponds to an attractive force per unit length ^f- 

Apart from reproducing this limit, Eq. (|16|) gives the explicit (exponentially 
decreasing for L >> W) corrections to this limit. 

Here, it is interesting to note that an alternative expression for the Casimir 
energy of nanotubes was given in Eq. (69) of Ref. 1261 Such expression, which is 
obtained through a combined use of the zeta function regularization and the gener- 
alized Chowla-Selberg formula, proves more useful in the limit W- oo. With our 



approach, such expression can be obtained by applying the inversion formula (JR 
to the sum over n (instead of the sum over I as we have just done) in Eq. (fT2"|) . Thus 
one gets, for the Casimir energy per unit length of the nanotube 



Ec 
W 



2.9 S 



LW 



(-iy 



6 Ann W 

r l 



9s 9v \ ^ cos (nw S) 



ttL 2 ^ 

n=l 



(18) 



This coincides with the result in Ref.[5Sl except for the fact that, in this reference, 
the spin degeneration hasn't been included, and there is a different sign (seemingly 
a misprint) in front of the double sum. 

In the limit ^ — )• oo (long nanotubes), one has = X^^Li cos ^ s ' . The 
corrections due to finite lenghts are also exponentially decreasing, except for the 
case of conducting nanotubes (5 — 0), where the I = term in the double sum is 
given by -^E. 



6. Final comments and conclusions 



Even though the experiments on graphene devices are still at a very preliminary 
stage, we have shown in this paper that MIT bag boundary conditions are the most 
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promising ones among the family of local boundary conditions that can reasonably 
(in the sense of our assumptions in section [2]) be imposed in the continuum limit. 
In particular, for nanoribbons, they predict not only the existence of a gap, which 
does not depend at all on the orientation, but also the existence of equally spaced 
energy levels. In the case of circular nanodots, they are the only ones, in the family 
of local boundary conditions studied here, predicting a spectrum which is symmet- 
ric around zero, as the experimenJIH seems to imply. Note that Ref. [55] had already 
suggested that this could be the case. In that reference, MIT bag boundary condi- 
tions are rather called Berry-Mondragon boundary conditions, because they were 
also studied in Ref. [19j in the context of quantum billiards. Indeed, this is the first 
guess a particle physicist would make when asked for confinement in a Dirac theory. 
Though experiments on graphene nanodevices are still at a very initial stage, the 
comparison with the experimental results presented in this paper shows this to be, 
most probably, the case. As for zigzag conditions, the fact that they fail to define a 
Fredholm Hamiltonian would be an obstacle to the formulation of a good quantum 
theory. 

In view of the previous considerations, we have performed in this paper a full zeta 
function calculation of the Casimir energy for nanotubes of arbitrary dimensions 
and chiralities, which lead us to an expression particularly adequate to the case of 
nanoribbons. Such expression clearly shows that, in this limit, and independently 
of the chirality of the initial nanotube, there is an attractive force between the 
boundaries of the ribbon, and that the corrections to this force decay exponentially 
with -pp. For some previous work on Casimir energies in the presence of graphene 
sheets, we refer the interested reader to Refs. [571 1551 
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